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1.  INTRODUCTION 


An  Important  quantity  obtainable  f»'om  Hatellltu  to  Hatolllto  tracking  ISfiT) 
iH  the  ] lne“Of -sight  accelerat l -d.  The  utilization  of  this  quantity  has  been 
attractive  to  geophysicists  especially  because  the  line  of -sight  acceleration 
can  be  expressed  almost  entirely  in  terms  of  the  gradient  of  the  disturbing 
potential  residual  to  a  chosen  reference  gravity  field.  In  previous  reports  and 
papers,  approximations  were  introduced  into  the  1 ine-of-slght  acceleration  model 
which  transformed  the  above  qualifier  "almost  entirely"  Into  "entirely",  If 
warranted,  such  approximations  are  of  great  practical  value  since  they  lend  to 
an  efficient  adjustment. 

Having  the  importance  of  a  rigorous  llne-of-slght  acceleration  model  in 
tnlnd,  we  address  the  issue  of  errors  in  the  app'-oximate  model  from  both  the- 
analytical  and  the  computational  standpoints,  the  latter  via  results  of  computer 
olfflulatlons.  The  emphasis  In  this  task  is  on  SST  In  the  high-low  mode,  n 
concept  that  has  grown  popular  among  geophysicists  as  Is  documented,  e,g,,  by 
Jekcll  and  Upadhyny  [1090],  These  authors,  as  well  as  Gleason  [1091],  envision 
the  "low"  satellite  as  the  space  shuttle  (or  a  satellite  at  the  shuttle's 
altitude),  and  the  "high"  satellite  as  one  (or  more)  of  the  OPS  satellites,  Due 
to  the  great  altitude  of  the  "high"  satellite,  the  gradient  of  the  disturbing 
potential  residual  to  a  sufficiently  detailed  reference  gravity  field,  such  as 
that  represented  by  a  (fi,6)  or  an  (8, ft)  spherical -harmonic  expansion.  Is 
negligibly  small  at  points  along  the  "high"  orbit.  This  allows  for  a 
simplification  whereby  the  gradient  of  the  disturbing  potential  may  be 
considered  only  at  observation  points  along  the  "low"  orbit.  However,  the 
derivations  in  the  present  study  will  proceed  without  this  simplification,  which 
will  he  introduced  only  an  a  last  step. 

The  ensuing  "last-ntep"  formula,  developed  later  on  as  equation  (31),  Is 
written  In  the  form 

dp  »  '0*'  <  . , ,  , 

Here  and  throughout,  the  notation  d  loetiMries  a  given  qiiantlfy  In  a  residual 
gravity  field,  l,e,,  the  difference  between  this  quantity  In  the  actual,  or 
complete,  gravity  field  (where  Its  symbol  has  no  superscript)  and  the  same 


((iiantlty  as  computed  In  a  reference  gravity  field  (where  lt»  jynbol  has  the 
jiiiperscrlpt  c).  Thus,  we  have 

f#  .10 

Sp  p  -  p  *  reeldual  llne^of-alght  acceJerat<'i  .. 

where  p  Is  the  llne-of-»lght  acceleration  In  the  complete  gravity  field  and  p® 
is  its  counterpart  In  the  reference  gravity  field;  and 

Mj  -  Xj  -  xj  . 

where,  in  the  Inertial  coordinate  system  (X,Y,Z),  * (X^ ,Y^ , Zj )  Is  the 

acceleration  vector  of  the  first  (low)  satellite  In  the  complete  gravity  field 
and  X^  is  its  counterpart  in  the  reference  gravity  field;  In  addition,  e°  Is 
the  1 ine-of-alght  unit  vector  as  computed  In  the  reference  gravity  field. 

In  following  the  derivation  presented  in  (Hajela,  1978],  Gleason  [1991] 
introduces  the  approximations  alluded  to  earlier,  and  obtains  the  result 
corresponding  to  the  above  “last-step"  formula  with  the  dots  removed  (cf, 
equation  3,  ibid. ) : 

(Jp  ^  . 

This  outcome  is  reached  also  by  Jekeli  and  Upadhyay  [1990],  whose  derivation 
agrees  with  that  In  [Rummel,  I960].  Thus,  all  four  references  just  mentioned 
arrive  at  the  same  approximate  model.  Gleason  [1991]  uses  it  in  the 
intermediate  form  (Implied  by  equations  3  and  10,  Ibid.): 

dp  -  -vJ[(X.V,Z)*e®  , 

where  ^(X,Y,Z)  la  the  gradient,  expressed  in  inertial  coordinates,  of  the 
disturbing  potential  at  the  location  of  the  "low"  satellite.  From  the  known 
time  t  associated  with  dp,  he  then  expresses  ^  in  earth-fixed  polar 
coordinates.  Upon  computing  a  priori  covariances  as  needed  for  the  collocation 
a<IJustfflcnt  model,  ho  uses  the  latter  to  obtain  an  optimal  solution  for  a  set  of 
<5^,,  the  mean  radial  components  of  the  surface  gravity  disturbances. 

Clearly,  the  accuracy  of  the  results  from  simulated  observations  as  In 
[Oleason,  1991],  and,  especially,  the  accuracy  of  the  gravity-field 
determination  where  the  observables  are  actual  1 Ino-of-slght  accelerations 
hinges  on  the  accuracy  of  the  mathematical  model.  In  order  to  design  a  rigorous 
procedure  for  adjusting  such  accelerations,  one  has  to  express  the  form  of,  and 


2 


(■valiiatw  the  effect  of.  the  terms  represented  by  the  dots  in  the  "last-step" 
formula  presented  at  the  outset.  Some,  perhaps  most,  of  these  terms  may  be 
fieKlijjihJe  In  general.  Others  may  be  negligible  upon  certain  restrictions  (on 
the  time  Interval  between  the  epoch  of  observation  and  the  epoch  at  which  the 
fjtate  vector  is  given,  on  the  minimal  degree  and  order  of  the  reference  gravity 
field,  etc,).  However,  if  there  exists  a  term  that  cannot  bo  made  negligibly 
small,  preferably  much  smaller  than  0.1  mgal  characterizing  the  magnitude  of  the 
noise  in  an  electronically-induced  signal,  this  term  should  be  a  part  of  the 
mathematical  model.  In  theory,  it  should  then  be  taken  into  account  In  the 
formation  of  observation  equations,  a  priori  covariances,  etc. 

Some  of  the  results  presented  by  Oleason  [1991]  Indicate  that  not  all  of 
the  terms  represented  by  the  dots  (discussed  above)  are  negligible.  When  he 
attributed  zero  errors  to  observables  simulated  in  one-second  intervals  within 
a  10*“squaro  equatorial  oceanic  area,  he  recovered  the  "ground  truth"  in  1* 
squares  with  no  better  than  a  1.1  mgal  r.m.s.  error  and  a  3.8  mgal  maximum  error 
(worse  results  wore  obtained  for  lower  observational  density).  In  performing 
comparable  simulations  over  a  larger  (iS*K20*)  and  geophysically  less  tranquil 
continental  area,  he  arrived  instead  at  a  1.0  mgal  r.m.s.  error  and  a  6,0  mgal 
maximum  error,  It  is  likely  that  the  cited  errore  are  Imputable  to  the  above 
modeling  approximations,  especially  upon  considering  the  magnification  of  errors 
caused  by  the  downward-continuation  problom.  Therein  lies  the  motivation  for  a 
(•('flnemcnt  of  the  traditional  1  Ine-of-slght  acceleration  model. 

The  present  study  concentrates  on  the  development  of  the  mathematical  model 
rigorous  to  within  the  first-order  differential  quantities.  Two  second-order 
dlfforeiitial  quantities  will  appear  (temporarily)  Incidental  to  the  derivations, 
while  further  second-  and  higher-order  differential  quantities  will  be  mentioned 
without  being  expressed  explicitly.  Although  the  "last-step"  formula  implies 
the  simplification  whereby  the  position  vector,  the  velocity  vector,  and  the 
acceleration  vector  of  the  "high"  satellite  in  the  complete  gravity  field  are 
cfjual  to  their  counterparts  computed  in  the  reference  gravity  field,  this 
.simplification  will  be  used  neither  in  the  derivations  nor  in  the  simulations. 
The  simulations  will  encompass  two  "low"  arcs  passing  over  the  above-mentioned 
oceanic  area,  each  limited  to  three  minutes  in  duration  from  the  epoch  of  a  i 
given  state  vector.  The  second  arc  will  be  used  in  conjunction  with  the 
reference  fields  represented  by  (6,0),  (8,8),  and  (10,10)  spherical -harmonic 


expnnslons,  whereas  for  the  first  arc  we  will  consider  only  the  (6,6)  reference 
field,  The  slwulatlons  will  serve  not  only  for  the  verification  of  the  derived 
mathematical  model,  but  also  for  the  assessment  and  possible  simplifications  of 
this  model  in  view  of  a  sufficiently  rigorous,  yet  computationally  manageable 
adjustment  of  the  line-of -sight  accelerations  in  the  high-low  mode. 

As  a  historical  note,  the  first  part  of  this  report  (Chapters  2-5)  was 
presented  in  a  similar  fora  in  a  GL  (APSC)  Internal  paper,  ”Satellite-to- 
Satfilllte  Tracking  in  the  High-Low  Mode:  Llne-of-Slght  Acceleration  In  a 
Roaidual  Gravity  Field",  April  1990,  Two  months  later,  a  OL  scientist,  Mr, 

David  M.  Gleason,  supplied  the  author  with  computer  simulations  supporting  the 
theoretical  findings.  An  expanded  version  of  the  paper  Including  a  summary  of 
the  computer  simulations  was  published  as  the  author's  Status  Report  No,  5  under 
the  present  contract,  period  4  April  1990  -  3  July  1990,  A  year  later,  another 
PL  (formerly  OL)  scientist,  Dr,  Christopher  Jekoli,  derived  an  equivalent 
algebraic  result  for  8p  by  taking  a  formal  differential  of  p.  However,  the  form 
of  this  result  is  not  considered  final,  A  similar  derivation  is  presented  in 
the  Appendix  herein,  where  the  differential  is  transformed  with  the  aid  of 
relations  from  Chapters  3  and  5  into  the  fin,il  formula  of  Chapter  5  (the  general 
version  before  the  last-step  simplification).  The  Appendix  is  instructive  in 
confl/'mlng  the  main  outcome  of  Chapter  5  via  a  semi -independent  route,  shorter 
overall  than  the  route  of  Chapters  2  5,  Other  points  of  discussion  are 

pro.seiiled  In  the  Appendix  Itaclf, 
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2.  GEOMETRIC  QUANTITIES  IN  THE  COMPEBTE  GRAVITY  FIELD 


The*  notation  In  t:hi»  analyala  1«  adopted  for  t ho  most  pari  from  [RummoJ , 
1960],  Thus,  the  inertial  axes  arc  denoted  by  the  letters  X,  V,  and  Z,  as  are 
the  natolllte  coordinates;  the  latter  arc  attributed  the  subscript  1  for  the 
first  satellite  (here  low),  and  the  subscript  2  for  the  second  satellite  (here 
high).  The  first-order  tine  derivative  is  identified  by  a  dot,  and  the  second- 
order  time  derivative,  by  a  double  dot.  Accordingly,  the  basic  symbolism  is 
presented  as  follows: 


.‘Jat,  1; 

2^1 

ft  (X,  , 

Zj), 

^1  * 

(X,  , 

Z,), 

- 

(X,  , 

V,,  z,): 

Sat,  2: 

h 

ft  (X2. 

'^2' 

Z2), 

rx2« 

^2- 

Zg). 

(X^, 

whore  the  underlined  quantities  are  vectors.  The  vector  differences  in 
position,  velocity,  and  acceleration  are  symbolized  by 


-12 


^^12 


^  2^1 


(1) 


As  will  become  apparent,  the  development  is  identical  whether  the  low-low  or  the 
high  low  configuration  Is  considered.  The  latter  case,  Ijowever,  allows  for  a 
almpl I f Icatlon  that  will  result  In  equations  {31)-(32d), 


The  quantities  that  may  be  subject  to  observation,  separately  or  in 
combination,  are  - 

p  ...  1 Inc-of 'Sight  range: 

I  e  dp/dt  1 ine-of-slght  velocity  (or  range  rate):  and 

...  2  2 

p  »  dp/dt  ^  d  p/dt  ...  line  of-sight  acceleration. 


All  of  the  above  symbolism  Implies  the  complete  (unabridged)  gravity  field. 

The  present  analysis  Is  concerned  with  p,  especially  In  view  of  expressing  the* 
latter  In  a  residua)  gravity  field.  The  1 Ine-of-slght  range  Is  given  by 

which  is  used  In  forming  e,  the  unit  vector  In  the  direction  Sat,  1  -  Sat,  2: 

1  V 

n  ^  h  X,2  ■ 


(3) 


The  notation  such  as  X  is  avoided  here;  Instead,  this  quantity  would  be 

2 

written  as  ^-12^  ^  similar  statement  applies  with  regard  to  any 

other  vector. 

The  time  derivative  of  p  follows  from  (2)  as 

^  hz%z  ■ 


<•  - 

Equation  (.3)  yields  the  time  derivative  of  e: 

i  .  (dp  Vdt)  5,^  .  p  ‘  ■  -p-"  P  x,2  .  p  ’  X, 


e  --  P  (Xi2  -  P  e)  .  (5) 

We  notice  that  e"e=0  so  that  the  vector  e  is  perpendicular  to  e.  From  (4) 
it  follows  that 

P  -  Xj2-e  .  .  (fia) 

which,  due  to  (5)  and  (4),  becomes 

P  --  Xj2’i  "  -  P^)  • 

As  a  matter  of  interest,  (6b)  can  be  given  in  an  alternate  form  In  writing 


-12  -12 


<^12-^>S  ^  ^12, n  ^  ^  Xj2.n  ’ 


where  component  of  X^^  perpendicular  to  e,  we  have 

-  S,2-X,2  ^  .  |X:,3  . 

and  (6b)  can  be  presented  as 


P  "  X  -e  ^  p  |X, 
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3.  GEOMETRIC  QUANTITIES  IN  A  RESIDUAL  GRAVITY  FIELD 


We  now  introducr?  a  reference  field  and  denote  the  quantities  computed  in  it 

c  •  r  ■  ■  c  c  *  c  '  (' 

by  a  superscript  c,  sucli  as  in  Xj  ,  X ^  ,  X^  ;  X^.  X^.  X ^  ;  and 


X^  -  x''  -  X*' 
-12  -2  -1  ■ 


•  P  *  P  •  P 


5?2  ■  ^2  -  S? 


In  analogy  to  (2)  and  (3),  one  has 


(X*"  •  x' 
^^12  ^12' 


-  ip"")  X'l^  .  (9) 

Precisely  the  same  procedure  that  resulted  in  (4),  but  carried  out  in  the 
reference  field,  yields 

•c  -c  c  , 

p  =  ^i2‘-  ■  ® 

For  the  "high"  satellite,  it  is  usual  to  adopt  and  X2  =  X2.  but 

here  this  simplification  will  be  left  for  the  very  end.  In  differencing  (4)  and 
( 9 ' ) ,  we  obta i n 

h  -  h""  -  Xj2-e  .  (10) 

The  difference  on  the  left-hand  side  of  (10),  which  owes  its  existence  to 
the  residual  gravity  field,  is  symbolized  by  6p: 

Sh  =  P  -  p‘^  ■  (11) 

This  quantity  Is  called  the  residual  lineof  sight  velocity.  In  order  to  bring 
(10)  to  a  convenient  form,  one  makes  the  following  approximation: 


Thus,  (10)  is  written  ,is 


where 


dp  ^j^.e 


^^12  ^12  ^-.2 


Although  the  symbol  -  has  been  I'ep laced  in  (13)  by 
In  (13)  should  be  kept  In  mind. 


the  approximation  present 


In  a  similar  vein,  the  procedure  that  resulted  in  (6b),  but  carried  out  in 
the  reference  field,  gives 

c  ^c  c  ,  c.-l,  ‘c  ,2  /•C.2, 

P  =  ^12*-  ^  ^'-12!  "  (P  )  1  •  05) 

In  subtracting  this  equation  from  (6b),  and  introducing  an  additional 
approximation  akin  to  (12),  namely 

P  ^  P^  .  (15') 


we  obtain 


where 


6p  =  "  (p'')‘^[|Xj2l^  -  IX 


12  ' 


-  P 


,  *0,2, 

(P  )  ] 


«5P  =  P 


c 

P 


^12  =  ^12  ■  ^?2 


(16) 


The  quantity  <5p  Is  the  residual  line-of-sight  acceleration.  Working  in  terms  of 
c  c 

e  ,  p  ,  etc.,  rather  than  in  terms  of  e,  p.  etc.,  has  the  practical  advantage  in 
that  the  reference  field  is  known  and  allows  for  the  actual  computation  of  such 
quantities.  We  re-write  (16)  as 


(16'  ) 


where  <5T  rcpresr*  the  quantity  inside  the  brackets  of  (16): 

•  2  •€  2  -2  *0  2 
<5T  |Xj2l  -  |Xj2l  -  P  ^  (P  ) 


(16"  ) 


At  a  later  stage,  the  approximations  (12)  and  (15')  will  be 
obtain  results  rigorous  to  the  first  order,  i.e.,  to  within 
differentials.  Corrections  due  to  both  approximations  will 
(16')  directly:  subsequently,  a  correction  due  to  (12)  will 


removed  in  order  to 
the  first 
affect  (16)  and 
affect  (5T  via  ( 13)  . 
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4.  REFORMULATION  oi  THE  RESIDUAL  LINE-OF-SIGHT  ACCELERATION 


We  consider  the  vector  ^^2  M4),  which  can  be  written  in  several 


forms : 


^12  =  ^12  ^12  =  ^2 


•  c  •  c 

(^2  -  > 


=  ^2  -  ^2 


(Xj  -  Xj  -  «2  -  • 


(17) 


where  the  ]ast  two  symbols  have  been  introdiired  as 

\c 


^1  • 


^2  ^  ^2  ?^2  • 


In  returning  to  the  formalism  of  (14).  we  have 
2 


1^12' 
and,  accordingly, 


^^12  ^^2>-'^,2  -  ^^2>  =  '^12''  2^12-^52  ‘  '4'" 


!X^2i 


1^12'' 


2^12-^12 


|X 


c  ,2 
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If  this  expression  is  substituted  into  (16"),  it  follows  that 


ft  ■ 


,’c  2 
(P  ) 


•  c  2  *0  •  c 

or,  upon  recalling  (14)  and  the  fact  that  IXjgl  ^^j2‘-22' 


(p 


c  ,2 


(18) 


We  first  denote  the  last  two  terms  in  (18)  by  6K ,  and  express  them  as 

aw  '2  ,  *c,2  ‘  ’C...'  ’C 

dK  =  -fl  +(p)  =  ( p  p  ) (p  +  p  )  . 

From  (11),  the  first  factor  on  the  right-hand  side  is  (5p  and  the  second  factor 

•  •  c 

is  <5p  +  2p  ,  yielding 

SK  -  <Sp^  2p^  Sp  . 

•  c  • 

Upon  adopting  p  from  (9')  and  i^p  from  (13),  which  contains  the  approximation 
alluded  to  earlier,  it  follows  that 

<5K  (5p^  2(X^2'5‘ 
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(  19  ) 


This  result,  is  only  intermediate,  since  It  does  not  account  for  the  effect  of 
-12^-12  -12  below.  Moreover,  it  contains  two  terms  (the  first 

two  terms  inside  the  last  parentheses)  corresponding  to  second  order 
differentials,  which  subsequently  will  be  neglected 


5.  RESULTS  RIGOROUS  TO  THE  FIRST  ORDER 

To  correct  two  out  of  throo  approximaf ions  implied  in  the  preceding 
development,  we  return  to  equation  (16)  and  observe  that  it  was  obtained  from 
(6b)  and  (15)  upon  adding  the  following  expression  to  (6b): 

V  ^  V  /Cv-l,*  ,2  *2,  i/iu  ,2  *2. 

-12*-  -12*-  ^  ^1-12^  P  )  -  p  (IXjgi  -  P  )  • 


Accordingly,  the  correction  (i.e.,  minus  the  above  expression)  is 


^  <5P  -  P^) 

(22) 

where  the 

first  term  is  due  to  (12), 

and  the  second  term  is  due  to  (15'). 

Upon 

consulting  (2),  one  forms 

^P  '  -12*-  ■ 

(23a) 

Sp'^  =  ^12*-  • 

(23b) 

With  the  aid  of  (23b),  the  differential  of  (3)  is  obtained  as 

where  (3)  has  also  served  in  a  substitution.  Since  the  second  term  inside  the 
brackets  is  the  vector  Projected  onto  e,  it  follows  that 

=  0  '  m,2,„  .  oe 

where  the  subscript  n  indicates  the  direction  normal  to  e.  We  recall  that 

0 

following  (19),  n  indicated  the  direction  normal  to  e  .  However,  this 

difference  in  interpretation  is  of  little  consequence  at  the  present  stage. 

In  particular,  the  quantities  e,  p,  etc.  in  (22)  (24)  can  be  replaced  at  will 
c  r 

by  e  ,  p\  etc  ,  since  such  approximations  introduce  errors  corresponding  only 
to  second  and  higher-order  differentials.  These  replacements  will  henceforth 
be  made  without  further  mention. 

Next,  with  the  aid  of  (4)  we  develop  the  last  factor  in  (22); 

Since  be  decomposed  into  two  orthogonal  constituents  in  a  familiar 


manner  as 


—12  -] 2  -  —  — 1 2 , n 


it  follows  that 


^12'^^12  -  ^^12-^»  *  ^12*^12. n  ’ 


and  (25a)  becomes 


IX, jl'  -  ■  j,2-S,2,„ 


(25b) 


The  first  term  In  (22).  which  owes  its  existence  to  (12),  is  denoted  c^ .  and  the 
second  term  in  (22).  which  is  due  to  (15'),  is  denoted  c^.  With  the  aid  of  (24) 


and  (25b)  together  with  (23b),  we  write 
^2  = 


12  — 12,n  ■ 


C.-2 


c  ,  ,  :,c  •  c 


C3  =  -(P  )  (^,2-')(Xj2-X,2.„) 


(26a) 

(26b) 


An  additional  correction,  denoted  c^,  is  due  to  the  substitution  of  6p 
from  (13)  into  <5K  in  (19),  which  entails  the  approximation  (12).  In  particular, 
in  obtaining  (19)  we  have  replaced  the  rigorous  value  of  6p  from  (10)  by  the 
expression  in  (13),  i.e.,  have  deformed  5p  by 

Xi2'e"  -  ^12-^ 

Accordingly,  the  correction  for  5p  will  be 
Xi2'^  • 

From  the  relation  preceding  (19)  it  is  apparent  that  the  correction  for  <5K  and 
thus  also  for  <5T  will  be 

and  from  (16')  it  follows  that  the  correction  for  5p  will  be 

-2(p'^)~^  ^12'  — 

In  consulting  (24).  we  finally  deduce  that 

■^4  -  -2(P')  ' 


1  2 


(27) 


We  are  now  In  a  position  to  formulate  the  final  version  of  6p,  rigoroua  to 
within  the  first  differentials.  In  returning  to  (21),  we  discard  the  terras 


which  are  essentially  second-order  quantities,  as  was  mentioned  earlier. 
(Consistent  with  this  fact,  Hajela  (1978],  Rummel  [1980],  and  Jekell  and 
Upadhyay  [1990]  show  that  the  above  expression  Is  negligible.)  Clearly,  all 
the  other  second-  and  higher-order  differential  quantities  generated  by  the 
replacements  described  following  equation  (24)  are  neglected  as  well.  Before 
transcribing  (21)  with  the  corrections  c^,  c^ ,  and  c^  Included,  we  present  the 
identity 


•  c  •  •  c  • 

-ax,,, 

-12  — 12, n  -12, n  — 12 


(28) 


which  shows  that  the  Index  n  can  be  shifted  from  one  constituent  to  the  other. 
This  fact,  easy  to  verify,  applies  In  conjunction  with  general  vectors,  ndt  only 
those  shown  In  (28);  here  It  will  be  used  also  for  c^  In  (26a)  and  c^  In  (27), 
What  remains  of  the  second  terra  In  (21)  will  be  denoted  as  Cj ,  and  will  be 
transcribed  as  Indicated  In  (28).  In  conjunction  with  c.  In  (26b),  we  will  use 
the  identity 


C  AC  _  AC 

12  -12,n  ■  -12.n  -12,n 


'4,nl 


which  again  holds  true  for  general  vectors. 


(28’) 


In  transcribing  and  completing  (21)  according  to  the  above,  we  have 


«P  -  gia-e  *  >  <^3  •  • 


(29) 


where  the  first-order  corrections  are 


<^1  -  4.n-ffl'2  ■ 

°2  -  i^?2.n-^12  • 

,  C.-2  .'C  ,2  c 

C3  =  (P  )  1^12. n^  -  *—12  ’ 

C4  =  -2(p  )  p  ■ 


(30a) 

(30b) 

(30c) 

(30d) 
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Since,  in  practice,  the  reference  field  is  sufficient  for  an  accurate 
description  of  the  position,  velocity,  and  acceleration  of  the  "high"  satellite, 
it  is  permissible  to  adopt 

^12  =  -^1  •  ^12  =  ’^1  •  ^12  "  “^1  • 

In  this  case,  the  formula  (29),  (30a-d)  simplifies  to  read 

Sp  =  -<5X  -e^  +  c;  ^  cl  +  cl  +  c’  (31) 

_ — I  ~ _ 1 _ 2 _ 3 _ 4 

where  the  first-order  corrections  are 


(32a) 

(32b) 

(32c) 

(32d) 
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6.  SIMULATED  RESIDUAL  LINE -OP-SIGHT  ACCELERATIONS. 
AND  ASSESSMENT  OF  THEIR  MODELING 


The  simulations  in  this  study  involve  two  arcs  of  the  "low”  satellite 
generated  by  Gleason  [1991].  The  first  arc  herein  is,  in  fact,  the  very  first 
of  143  generated  arcs  entering  the  10“-square  equatorial  oceanic  area  bordered 
by  the  lines  of  5°  E,  W,  N,  and  S.  And  the  second  arc  herein  is  the  arc 
associated  with  the  greatest  modeling  errors  among  the  143  arcs;  below  we  refer 
to  it  as  the  "worst"  arc.  The  "low"  satellite,  denoted  Sat.  1.  emulates  a 
typical  space  shuttle  orbit  of  approximately  300  km  altitude,  while  the  "high" 
satellite,  denoted  Sat.  2,  emulates  one  of  the  GPS  satellites  of  approximately 
20,189  km  altitude.  Among  the  GPS  satellites,  the  one  chosen  for  simulations  of 
the  residual  line  of-slght  accelerations  is  the  one  which  has  a  minimal  zenith 
distance  with  respect  to  the  "low"  satellite.  In  [Gleason.  1991]  the  average  of 

such  minimum  distances  is  reported  to  be  28°.  Although  the  duration  of  the 

"low"  arcs  over  the  oceanic  region  of  interest  is  two  minutes  or  less,  the  two 
arcs  serving  in  the  present  study  are  extended  to  181  s.  Gleason  [1991]  uses 
Rapp's  1981  geopotential  coefficients  through  harmonic  degree  and  order 
(180,180)  to  describe  a  field,  which  he  calls  "true",  approximating  the  complete 
gravity  field  of  the  earth,  while  for  the  reference  gravity  field  he  chooses  a 
(6,6)  subset  of  these  coefficients. 

The  mathematical  model  for  the  residual  line  of-sight  acceleration  is  (29), 
(30a-d).  The  above  first  arc  serves  mainly  for  the  determination  of  the  step 
size  used  in  the  Hamming  orbit  generator,  which  would  yield  the  time  derivatives 

with  sufficient  accuracy  for  any  epoch  within  three  minutes  from  the  epoch  of 

the  state-vector.  Although  in  [Gleason.  1991]  the  step  size  is  one  second,  here 
three  step  sizes  have  been  tested:  At  1  s.  At  0  1  s.  and  At^O  01  s  At  a 
desired  epoch,  the  time-derivative  is  computed  by 

[value (epoch ^  At  )  -  value (epoch  At )]/ 2At 

The  step  size  At  1  s  is  deemed  too  coarse  for  the  precise  calculations 
needed  in  the  model  verification.  In  particular,  errors  in  6p  have  been  found 
to  range  from  0.04  mgal  to  0.18  mgal  for  epochs  ranging  from  2  s  to  179  s  from 
the  state  vector  epoch.  The  step  size  At  =  0 . 1  s  reduces  the  magnitude  of  these 


errors  by  over  90*!.  To  produce  results  accurate  to  0,001  mgal  or  better, 

At “0,01  8  Is  adopted.  The  step-size  errors  affect  the  "true"  dp  and  the  first 
term  on  the  right-hand  side  of  (29),  called  "basic  term",  by  approximately  the 
same  amount.  Thus,  the  error  of  the  basic-term  modeling  (l.c,,  the  difference 
between  the  basic  term  and  dp),  as  well  as  the  error  of  a  more  complete 
modeling,  are  affected  by  the  step-size  errors  to  a  much  lesser  degree  than 
0.001  mgal,  perhaps  by  two  or  more  orders  of  magnitude.  We  conclude  that  for 
the  modeling-error  analysis,  the  step  size  dt-0.01  s  Is  entirely  satisfactory. 

With  this  step  size  and  a  (6,6)  reference  field,  the  values  of  dp  on  the 
first  arc  have  been  generated  in  four  Intervals  from  the  state-vector  epoch: 
28,  60  8,  120  8,  and  179  s.  These  four  epochs  will  be  used  throughout  the 
analysis.  The  simulated  values  of  dp  at  these  epochs  are 

-2.883  ,  -3,370  ,  -1.697  ,  3.782  , 

the  units  being  mllllgals.  The  corresponding  values  of  the  basic  term  differ 
from  the  dp  values  by  the  amounts,  called  errors,  which  are  listed  (In  mgal)  as 

0,002  .  0.016  ,  -0.043  ,  -0,142  .  (33a) 

When  the  term  Cj  is  Included,  the  errors  decrease  substantially: 

■  ,  -0,0004  ,  -0,0009  ,  0,0011  ,  (33b) 

Finally,  when  also  the  terms  c^.  c^,  and  arc  Included,  the  errors  become 

--  ,  -0.000005  ,  0,000003  ,  0,000002  ,  (33c) 

The  notation  " — "  Indicates  completely  negligible  amounts.  Already  at  this 
early  stage,  the  validity  of  the  mathematical  model  (29),  (OOa-d)  is  confirmed. 
The  small  remaining  errors  displayed  In  (33c)  stem  partly  from  neglecting  the 
second-  and  higher-order  differentials,  and  partly  from  the  step-size  errors: 
tiie  round-off  errors  are  negligible  here  (with  the  double-precision  arithmetic) 

At  a  state-vector  epoch  the  errors  are  zero,  since  the  position  vector 
and  the  velocity  vector  in  the  complete  gravity  field  are  equal  to  their 
counterparts  in  the  reference  field  by  construction.  (For  the  "high"  satellite 
this  equality  holds  true,  or  nearly  true,  also  at  epochs  different  from  the 
state-vector  epoch.)  Accordingly,  at  a  state-vector  epoch  we  have 


which  implies  that 

^12  “  '  ^12  “  0  = 
e  p  .»  . 

Thus,'  all  four  corrective  terms  c,  through  c.  Jn  (29),  as  well  as  all  second- 
and  higher-order  differential  terms,  are  zero,  and  the  basic  term  supplies  the 
rigorous  value  at  any  epoch  where  the  state-vector  is  given.  The  main  question 
to  bo  answered  is:  How  far  from  the  state-vector  epoch  is  the  basic  term 
accurate  enough  to  represent  alone  an  acceptable  mathematical  model? 

In  considering  the  first  arc,  and  requiring  the  modeling  accuracy  to  be  at 
least  0.1  mgal,  we  observe  from  (33a)  the  basic  terms  alone  would  be  acceptable 
for  epochs  separated  by  nearly  three  minutes  from  the  state-vector  epoch, 
However,  the  situation  depicted  by  the  first  arc  is  overly  optimistic.  If  we 
consider  instead  the  “worst"  arc  in  the  same  context,  l,c,,  with  a  0,0i  s  step 
Hl/,e  and  a  (6,6)  reference  field,  in  lieu  of  (33a)  we  obtain 

0,010  ,  0,269  ,  0,386  ,  0,545  ,  (34) 

Bince  the  pertinent  oceanic  area  is  geophysically  tranquil,  the  relatively  large 
errors  in  (34)  may  be  the  norm  for  many  other  regions.  Equation  (34)  Indicates 
that  the  basic  term  alone  could  represent  a  (marginally)  acceptable  model  only 
in  the  case  of  observational  epochs  closer  to  the  state-vector  epoch  than  30  s. 

We  now  present  in  detail  the  results  that  culminated  in  (34).  This  will 
confirm,  once  again,  the  validity  of  the  model  (29),  (30a-d),  and  suggest  a 
simplified  residual  1 ine -of-sight  model,  although  not  as  simple  as  the  basic 
term  alone.  We  recall  that  the  step  size  is  0,01  s,  the  "true"  gravity  field 
is  given  by  (180,160)  spherical-harmonic  coefficients,  the  reference  field  is 
given  by  the  (6,6)  subset  of  these  coefficient,  end  the  results  are  listed  in 
milligals.  Two  additional  groups  of  results  will  be  presented,  corresponding 
to  iB,8)  and  (10,10)  reference  gravity  fields.  As  peripheral  information,  we 
list  the  values  of  8p  for  the  throe  reference  fields  as  obtained  for  the  usual 
four  epochs  (2  s,  60  s,  120  s,  and  179  s  from  the  state-vector  epoch): 
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(6,6)  2.607  ,  -2.000  .  0.417  ,  2.038  ; 

(8,8)  3.077  ,  -0.423  ,  3.314  ,  6.191  ; 

(10,10)  .  -0.530  ,  -5.069  ,  -1.788  ,  1.130  . 


These  values  do  not  reveal  any  clearcut  pattern,  or  significant  differences  in 
magnitude,  which  would  single  out  one  reference  field  in  preference  to  another. 

In  proceeding  in  detail  with  the  (6,6)  reference  gravity  field,  we  will 
list  all  the  results  to  six  decimals.  The  term  c^  will  be  presented  separately, 
while  the  terms  c  c  ,  and  c  which  all  depend  on  M  ,  will  be  combined  into 

^  o  4  i  fa 

the  term  *^2  *  ^3  ^  ^^4  sake  of  Interest,  a  term  containing  the  second- 

order  differential  quantities  seen  in  (21)  will  also  be  listed,  written  as  (z); 

(z)  =  (p^^)  ■ 

Clearly,  this  term  accounts  only  for  some  second-  and  higher-order  differential 
quantities.  Its  listing  will  illustrate  that  it  is  negligibly  small,  as  already 
transpired  from  fflajela,  I978J.  [Rummel ,  1980],  and  [Jekeli  and  Upadhyay,  1990]. 

The  guidance  to  the  foui'  groups  of  results  associated  with  the  "worst"  arc 
is  as  follows  Above  each  group  is  listed  the  time  interval  from  the  state- 
vector  epoch  (2  s,  60  s,  120  s,  or  179  s),  together  with  the  corresponding  value 
of  <5p  (see  also  the  preceding  outcome  for  the  6.6  reference  field,  where  these 
values  are  shown  to  three  decimals).  Each  group  proper  begins  with  the  line 
headed  "individual  terms",  where  the  first  item  is  the  basic  term,  the  second 
item  is  the  term  c^,  the  third  item  is  the  combined  term  ^2^^3^^4’ 
fourth  item  is  (z)  The  next  line  is  headed  "modeled  <5p" .  where  the  model  is 
presented  in  a  cumulative  fashion;  the  first  item  is  the  repetition  of  the  basic 
term,  the  second  item  is  [basic  term  +  c^].  and  the  third  and  last  item  is 
[basic  term  ^  ^  ^3  ^ ^4  ^  ^ complete  right-hand  side  of  (29). 

(The  value  in  the  second  or-  the  third  item  is  found  as  the  value  in  the 
preceding  item  on  the  same  line  plus  the  value  just  above  the  desired  item.) 
Finally,  the  third  line  is  headed  "remaining  errors",  and  it  contains  the  errors 
in  the  values  representing  "modeled  dp"  of  the  second  line.  (Each  of  the  three 
items  on  the  third  line  is  computed  as  the  value  Just  above  it  minus  the  value 
<5p  listed  above  the  first  line  )  All  values  are  rounded-off  to  six  decimals. 
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The  four  groups  of  results  are  listed  below. 


2  s 

.  Sp  -  2. 

607035 

basic  term 

'^1 

(z) 

individual  terms: 

2.616992 

0.009953 

0.000004 

(  --  ) 

modeled  6p: 

2.616992 

2.607039 

2 . 607036 

remaining  errors: 

0.009957 

0.000005 

0.000001 

(35a) 

60  s 

.  =  -2 

.000465 

basic  term 

‘^2*‘'‘3^"4 

(z) 

individual  terms: 

-1 .731951 

0.265731 

0 . 002781 

(0.000001 ) 

modeled  5p: 

1 .731951 

1 .997682 

-2 . 000463 

rema i n i ng  errors : 

0.268514 

0.002783 

0.000002 

(35b) 

120  s 

.  <5p  =  0. 

416935 

basic  term 

^1 

(z) 

individual  terms. 

0.803185 

-0.378171 

-0.008090 

(0.000004) 

modeled  6p: 

0  803185 

0.425013 

0.416924 

remaining  errors: 

0.386250 

0.008078 

0 . 000012 

(35c) 

179  s 

.  <5p  =  2 

037552 

basic  term 

^2*‘^3*^4 

(z) 

individual  terms: 

2.582464 

0  530705 

0.014204 

( -0.000001 ) 

modeled  6p\ 

2.582464 

2  051759 

2 . 037555 

remaining  errors: 

0.544912 

0.014207 

0 . 000003 

(35d) 

1  9 

For  the  sake  of  interest,  we  list  the  six  state-vector  components  (in 
inertia]  coordinates)  at  the  epoch  179  s  for  the  "low"  satellite,  followed  by 
the  six  state  vector  components  at  the  same  epoch  for  the  "high"  satellite,  both 
state-vectors  having  been  generated  in  the  "true"  field;  the  position  components 


are  given  in  meters,  and  the 
Sat.  1:  X  695,616  623  . 

X  -  -6,772.261608 

Sat .  2:  X  -  10,970,235  262 

X  =  -1,651. 398895 


velocity  components  are 

Y  6,628,710.831  . 

Y  612  098398  , 

Y  =  22,255,878.282  , 

Y  =  2,029.425766  . 


given  In  meters  per  second. 

Z  -  -190,943.856  , 

Z  =  -3,683.667565  ; 

Z  =  -9,466,691.250  , 

Z  =  2,857.380367  . 


Accordingly,  at  this  epoch  we  have  1X^1=6,667,844  m,  i.e.,  the  altitude  of  the 
first  satellite  is  about  290  km;  and  1X21=26,557.267  m,  i.e.,  the  altitude  of 
the  second  satellite  is  about  20,180  km.  For  the  first  satellite,  the  component 
differences  and  are 

-0.081  ,  -0.450  ,  0.743  ;  0.000370  ,  -0.004181  ,  -0.007800  , 

The  first  three  values  reveal  that  1=0.872  m.  For  the  second  satellite  the 

component  differences  are  essentially  zero  (they  are  0.000001  and  -0.000001  in 
X  and  Y.  respectively,  smaller  otherwise).  This  is  consistent,  for  example, 
with  the  statement  that  followed  equation  (30d). 

The  error  values  in  (35a-d)  warrant  a  discussion.  We  begin  with  the  third 
item  in  each  of  the  lines  (35a-d),  i.e.,  with  the  error  remaining  in  the 
expression  [basic  term  ^  ^  *"2  * ''O  *"4  ^  ^  error  is  exceedingly  small  for 

all  four  epochs,  confirming  the  theoretical  formula  and  indicating  that  the 
second  and  higher-order  differentials  are  negligible  within  three  minutes  from 
the  state-vector  epoch  (and  probably  much  beyond  three  minutes),  that  the  step- 
size  errors  are  inconsequential,  and  that  the  orbital  simulations  are  flawless. 
In  fact,  the  reliability  and  accuracy  of  the  orbit  generating  program  perfected 
by  Gleason  [1991]  is  verified  here  in  an  independent  and  interactive  fashion 
(interactive  in  the  setrse  of  high-low  satellites,  high-low  order  gravity 
fields).  In  continuing  with  the  second  item,  i.e.,  with  the  error  remaining  in 
the  expression  [basic  term  ♦  c^ ] ,  we  note  that  this  error  is  at  or  below  the 
level  of  0,01  mgal  for  all  four  epochs.  Thus,  the  corrective  term  c^  is 
sufficient  to  extend  the  validity  of  the  basic  term  to  three  minutes  from  the 


state-vector  epoch,  and  much  beyotut  three  minutes  if  one  is  willing  to  accept 
modeling  errors  on  the  order  of  0.1  mgal.  Finally,  the  first  item,  i.e.,  the 
error  in  the  basic  term,  indicates  that  the  model  comprising  the  basic  term 
alone  introduces  inadmissible  errors  already  within  a  few  tens  of  seconds  from 
the  state  vei^tor  epoch,  even  If  one  tolerates  the  modeling  error  of  0.1  mgal 

Faced  with  the  unacceptabi 1 1 ty  of  the  basic  term  for  accurate  modeling  of 
the  1 ine-of-slght  accelerations  in  the  gravity  field  residual  to  a  (6,6) 
reference  field,  one  is  compelled  to  ask  whether  a  reference  field  higher  than 
(6,6)  could  at  least  partly  remedy  this  shortcoming.  Accordingly,  two 
additional  computer  runs  have  been  analyzed,  utilizing  an  (“.8)  and  a  (10,10) 
reference  field.  First  of  all,  we  remark  that  the  errors  in  the  model  [basic 
term  ^  c^  *  ^  *^'2  *  ^  ^'4  ^  ^  again  exceedingly  small,  listed  as 

(8.8)  .  0.000005  ,  0.000001  ,  ,  0.000007  . 

(10,10)  ....  0.000006  .  0.000003  ,  0.000002,  0.000005  . 

There  is  no  essential  difference  between  these  errors  and  the  errors  presented 
above  for  a  (6.6)  reference  gravity  field. 

Next,  we  present  the  errors  of  the  basic  term  in  conjunction  with  these 
reference  fields,  again  for  the  four  chosen  epochs  (2  s,  60  s,  120  s,  179  s) 

We  enhance  the  usefulness  of  the  results  by  adding  also  the  errors  of  the  model 
[basic  term  »  parentheses.  These  two  kinds  of  errors  are  listed  below 

to  three  decimals.  For  the  sake  of  completeness,  the  errors  presented  erirlier 
for  a  (6,6)  reference  field  (see  equations  35a-d)  are  recapitulated  using  the 
same  format : 

(6,6)  .  O.CIO  (  )  ,  0.269  (0.003)  .  0.386  ( 0  008 )  .  0  (0  014)  : 

(8,8)  .  0,005  (--)  .  0.108  (0.002)  .  0.046  (0.005)  ,  0.041  (0  010)  ; 

(10,10)  ....  0,008  (  )  ,  0.188  (0.001)  ,  0.166  (0.001)  ,  0  162  (  ) 

Although  we  notice  improvements  in  the  basic  term  in  conjunction  with  higher 

order  referinice  fields  (as  compareti  to  the  6,6  field),  the  improvements  are  not 
monotonous.  For  example,  the  error  of  0.386  mgal  for  the  (6,6)  reference  fiehi 
improv<>s  to  0.046  mgal  for  the  (8.8)  reference  field,  but  then  it  grows  to  0  166 
mgal  for  the  (10,10)  reference  field. 


From  the  above  results  it  appears  that  a  re'atlvely  high-order  reference 
gravity  field  would  have  to  be  utilized  to  substantially  lower  the  error  level 
of  the  basic  term,  if  this  is  indeed  possible.  As  another  attempt  to  make  the 
errors  in  the  basic  term  acceptable,  one  could  lower  the  time  interval  between 
the  state-vector  epoch  and  any  observational  epoch  to.  eg.,  less  than  30  s. 
(Although  in  the  case  of  the  8.8  reference  field  the  errors  in  the  basic  term 
are  seen  to  decrease  beyond  the  epoch  60  s.  this  feature  is  exceptional  as  is 
gleaned  from  the  results  for  the  10.10  and.  especially  the  6.6  reference  field.) 
Both  these  possibilities,  whether  exploited  separately  or  in  combination,  would 
limit  the  practical  usefulness  of  the  mathematical  model  represented  by  the 
basic  term  alone.  Moreover,  such  a  model  would  still  suffer  from  poor  or 
marginal  accuracy. 

We  conclude  the  simulation  analysis  by  mentioning  that  due  to  computer 
limitations,  one  is  naturally  concerned  about  the  possibility  that  the  simulated 
errors  may  be  overly  optimistic.  This  possibility  stems  from  the  fact  that  the 
"true"  gravity  field  has  been  associated  here  with  a  (180,180)  field  rather  than 
with  a  (360.360)  or  a  more  detailed  field.  However,  we  alleviate  such  concerns 
by  presenting  the  errors  corresponding  to  a  (of)  reference  field  (as  used 
earlier)  and  to  a  mere  (60,60)  "tr  field: 

0.010  (--)  ,  0.26S  (0.003)  ,  0.380  (0,008)  ,  0.538  (0.014)  . 

While  (5p  at  any  of  the  four  epochs  has  changed  by  0.2  mgal  or  less  between  the 
(60,60)  and  (180.180)  "true"  fields,  the  errors  in  the  basic  term  are  seen  to 
have  changed  by  0.007  mgal  or  less,  and  the  errors  in  [basic  term  +  c^ ]  are  seen 
to  be  unchanged.  By  far,  most  of  the  change  in  the  residual  1 ine-of-sight 
acceleration  due  to  the  change  in  the  "true"  field  is  absorbed  by  the  basic 
term.  These  considerations,  extrapolated  to  moie  detailed  "true”  fields, 
indicate  that  the  simulated  error  analysis  is  lik(!ly  to  be  valid  with  sufficient 
accuracy  also  in  the  actual  gravity  field  of  the  e.irth 


7.  RECOMMENDATIONS  FOR  RESIDUAL  LINE-OP-SIGHT  ACCELERATTOfj  MODELING 


Guided  by  tlie  above  results  aud  considerations,  we  envision  a  tnanaKeable. 
yet  accurate  mathematical  model  lepresented  by  (basic  term  +  c^].  We  have  seen 
that  to  within  three  minutes  from  the  state-vector  epoch,  and  probably  much 
beyond  three  minutes,  the  errors  in  such  a  model  are  negligible  even  for  a  (6,6) 
reference  gravity  field.  Two  avenues  open  quite  naturally  before  us.  The  first 
is  to  use  the  term  c^  as  a  f.illv  equivalent  partner  of  the  basic  term  in  the 
model,  involved  in  the  formation  of  the  design  matrix  for  the  collocation 
adjustment,  of  a  priori  covariances,  etc. 

The  second  avenue  is  much  simpler  and  much  more  efficient.  In  consulting 

the  detailed  results  listed  in  conjunction  with  the  (6,6)  reference  gravity 

field,  we  observe  that  for  the  most  part  the  term  c^  is  significantly  smaller 
than  the  basic  term.  This  is  all  the  more  so  for  shorter  time  intervals  from 
the  state-vector  epoch,  and  for  higher  order  reference  fields.  Thus,  under 
certain  stipulations  (which  should  be  subject  to  furthei  ests,  including  tests 
covering  additional  regions  of  the  globe),  the  term  c^  could  be  treated  not  as 
a  full  partner  of  the  basic  term,  but  merely  as  a  small  observational  correction 
separate  from  the  model.  This  conceptual  simplification  is  far-reachiiig .  It 
makes  it  possible  to  revert  to  the  basic  term  as  a  highly  efficient  mathematical 
model  for  the  residual  1 ine-of-sight  accelerations  by  overcoming  the  initial 
shortcoming  of  this  model,  the  low  accuracy. 

The  crucial  aspect  of  this  approach  is  the  computation  of  the  correction 
Cj .  In  practice,  c^  can  be  approximated  by  the  same  procedure  that  has  been 

used  in  this  analysis.  The  true  gravity  field  can  be  approximated,  eg.,  by  a 

(180,180)  field,  and  the  satellite  orbits  in  question  can  be  approximated  by 
simulated  orbits  sufficiently  close  to  their  true  counterparts.  When  necessary, 
the  approximating  orbits  can  be  rectified  (to  correct  for  their  drift,  in  time, 
with  respect  to  the  true  orbits).  The  step  size  can  be  lowered  to  0.1  s,  etc 
(for  better  accuracy),  but  only  in  the  immediate  vicinity  of  the  obsei '-at iona 1 
epoch;  elsewheie  the  step  size  can  be  much  coarser  in  analogy  to  the  technique 
utilized  by  Gleason  [1991]. 

Akin  to  a  previous  concern  about  the  validity  of  error  analysis  using  a 
truncated  representation  of  the  gravity  field,  one  may  now  be  concerned  about 
the  accuracy  of  c^  computed  in  such  a  representation.  It  is  assumed  that  the 
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spherical-harmonic  coefficients  providing  this  representation  are  reasonably 
accurate.  An  example  using  a  (60,60)  set  of  coefficients  versus  a  (180,180) 
set,  where  the  former  Is  a  subset  of  the  latter,  again  helps  to  alleviate  such 
concerns,  in  writing  Cj  computed  with  a  (60,60)  truncated  field  In  brackets, 
and  -tn  writing  c^  computed  with  a  (180,180)  "control"  field  In  braces,  we  have 
for  the  usual  four  epochs: 

[-0.010]{-0.010},  [-0.266]{-0.266}.  [-0 . 372] {-0 . 378} .  [-0 . 524] {-0.531 } . 

Accordingly,  for  an  observational  epoch  separated  by  three  minutes  from  the 
state-vector  epoch,  the  effect  of  the  Indicated  truncation  generates  a  mere 
0.007  mgal  difference  In  the  value  of  c^,  corresponding  to  the  outcome  seen  at 
the  close  of  the  preceding  section.  We  reiterate  that  by  far,  most  of  the 
change  In  the  residual  1 Ine-of-slght  acceleration  Is  absorbed  by  the  basic  term. 
This  behavior  of  the  basic  term  further  supports  the  concept  of  using  c^  as  an 
observational  correction  computed  In  a  reasonably  detailed  and  reasonably 
accurate  "true"  gravity  field. 

In  the  above  setup,  an  approximation  to  the  true  correction  c^  can  be 
obtained  according  to  (30a)  or,  in  practice,  according  to  (32a).  In  principle, 
the  suggested  approach  uses  an  "adjoint  model"  (operating  on  assumed  satellite 
orbits  in  conjunction  with  an  assumed  gravity  field)  alongside  the  "true  model" 
(operating  on  the  actual  satellite  orbits  governed  by  the  actual  gravity  field), 
but  only  insofar  as  the  correction  Cj  is  concerned.  The  advantage  of  this 
approach,  based  on  the  In-practlce  ascertainable  smallness  of  c^  relative  to  the 
basic  term,  Is  substantial.  Due  to  Cj  being  considered  as  an  observational 
correction,  the  need  for  its  treatment  In  view  of  the  design  matrix,  of  a  priori 
covariances,  etc..  Is  eliminated. 
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8.  CONCLUSION 


In  the  first  part  of  this  study,  the  residual  1 i ne-of-slght  acceleration 

for  a  general  satellite  configuration  la  developed  rigorous  to  within  the  first- 

order  differentials.  In  addition  to  the  usual  "basic  term"  treated  in 

geophysical  literature,  two  kinds  of  corrective  terms  are  derived;  the  term  c^ 

(related  to  satellite  velocities)  and  the  combined  term  c„+c_+c.  (related  to 

2  3  4 

satellite  positions).  Subsequently,  the  general  formulation  is  specialized  for 
a  high -low  configuration.  In  this  case,  the  term  c^  depends  essentially  on  the 
residual  velocity  vector  of  the  "low"  satellite  at  the  epoch  of  observation,  and 
the  term  (c.+c  -►c  )  depends  essentially  on  the  residual  position  vector  of  the 

£  o  4 

same  satellite  at  the  same  epoch.  The  above  qualifier  "residual"  Implies  the 
difference  between  a  vector  expressed  in  the  complete  gravity  field  and  its 
counterpart  obtained  in  the  reference  gravity  field. 

In  the  latter  part  of  the  study,  the  analysis  of  computer  simulations  for 
a  high-low  configuration  confirms  the  first-order  general  formula  (devoid  of  - 
approximations  involving  the  "high"  satellite).  In  all  cases  examined  (two  arcs 
of  up  to  three  minutes  in  duration,  three  different  reference  gravity  fields), 
the  first-order  results  agree  with  the  known  "true"  values,  obtained  with  the 
step  size  At=0.01  s,  to  within  0.00001  mgal .  The  analysis  clearly  Indicates 
that  the  basic  term  lacks  the  accuracy  to  represent  alone  a  valid  mathematical 
model:  in  one  realistic  simulation  its  error  surpasses  0.5  mgal.  However,  when 
this  term  is  corrected  by  Cj,  its  accuracy  improves  significantly,  attaining  the 
level  or  0.01  mgal  or  better.  Thus,  the  term  (c^+c^+c^)  as  well  as  all  of  the 
higher-order  differentials  can  safely  be  neglected  in  practice.  As  a  further 
possible  simplification,  it  is  recommended  that  the  term  c^  be  considered  merely 
as  an  observational  correction  in  the  subsequent  (collocation)  adjustment, 
computed  in  an  "adjoint  model"  where  the  true  gravity  field  and  the  true 
satellite  orbits  are  substituted  for  by  their  sufficiently  accurate  and  easily 
manageable  approximations. 
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APPENDIX 


ALGEBRAIC  DERIVATION  AND  A  USEFUL  TRANSFORMATION 
OF  THE  FORMULA  GIVING  THE  RESIDUAL  LINE-OF-SIGHT  ACCELERATION 

To  find  the  formal  differential  of  the  line-of  sight  acceleration,  we  first 
recapitulate  the  standard  relations  (2).  (3),  (4),  and  (6b): 

P  ■  15,21  .  e  .  p-‘  5,2  •  P  ■  S,2-e  . 

P  ■  5,2-E  •  P‘'<!S,2l^  -  p'^ 

<2  -2 

In  forming  the  differential  Sp  of  p.  use  is  made  of  6p,  ,  and  6p  ; 

ip  -  OIS, 2-5,2)'''  ■  .  e.4X,2  ■ 

XIX, 2I'  ■  X(X„,.X,2)  ^  25,2-M,2  . 

6'p^  -  2p  Sp'  =  2‘p(Sk^^-e  >  Xj2-^) 


We  now  readily  deduce 


<5^  ^ 

X5,2-P-  •  5,2-Xo 

-  p-^  e.M,2 

(IX, 21"  -  p') 

where 

•  "  '‘^5,2-XX,2 

2p{SX^^-e 

^  Xj2-^)]  . 

-2 

- 1  - 1 

-2 

de  = 

P  <5P  Xj2  ♦  P 

^-^12  =  ^ 

^Xj2  -  P  (e-«, 2)^12  • 

Since  the 

differential  Sp 

pertains  to 

the  reference  field  identified  by  the 

superscript  c,  all  finite  quantities  above  are  attributed  this  superscript: 

where 

The  formulation  (A. 1,2)  is  not  very  Insightful  or  easily  tractable.  Moreover, 

■■  c 

we  wish  to  present  "p  in  a  form  where  the  "basic  term"  followed  by 


the  most  important  of  the  remaining  terms.  The  analysis  of  Chapter  6  indicates 
that  such  a  term  contains  ^  . 


In  the  general  formula  (29),  (30a-d),  this  term  is  presented  coiu  is(>ly  as 
.  Thus,  we  attempt  to  transform  (A. 1,2)  Into  a  form  such  as  (29),  (30a  d) 
First,  we  use  equation  (9');  tlie  decomposition  of  a  vector  into  two  orthogonal 
constituents  as  in  the  step  leading  to  equation  (24);  and  the  identity  (28): 


^f2-^12 


•  c  c  .  • 
p  e  »<5X 


12 


Xj2-f^l2  - 


^12-^12, u  ^  ^12,n-^12 


As  a  result,  we  can  already  e,xtract  the  term  c^  from  (A.l): 

Cl  -  2(p^)  ^  x;2,n-^12  • 


(  A  .  3  ) 


After  this  extraction,  only  the  third  term  remains  inside  the  second  pair 
of  brackets  iti  (A.l);  it  contains  which  is  also  featured  by  the  term 
following  the  basic  term  in  (A.l).  Before  proceeding  further,  we  give  ^  a  form 
more  convenient  than  (A. 2).  In  utilizing  (9)  in  (A. 2),  we  have 

The  expression  inside  the  brackets  is  that  referenced  in  the  preceding  paragraph 
with  regard  to  the  decomposition  of  a  vector  It  thus  follows  immediately  that 

This  relation  could  have  been  transcribed  directly  from  (24). 


Next,  we  express  the  second  term  on  the  right  hand  side  of  (A.l) 
(A. 4)  as  well  as  the  identity  (28)  and  the  three  lines  following  it 


We  use 


Xj2-^  ^  ^  ^?2-^12.n 


'  4,n-^12 


Thus,  the  desired  term  is 


’  ^^2.n-^12  • 


(  A  .  .S  ) 


To  express  the  third  term  on  the  right  hand  side  of  (A.l),  we  use  equation 
(9');  th(?  notion  of  a  vector  decomposition  utilized  already  twice  above,  and 
equat ion  ( 28 '  )  : 


•  r  •  c  c  c, 

-  <^12 ^ 

•  f'  •  p  •  p  2 

^12<2,n  =  1^12. n'  ' 


Thus,  the  desired  term  is 


,  C.-2  -c  ,2  c 

(p  )  |Xi2  n!  i  '^12  • 


(A, 6) 


The  remaining  term  (besides  the  basic  term)  on  the  right-hand  side  of  (A.l) 
is  that  mentioned  earlier,  namely  the  last  term  inside  the  second  brackets.  To 
give  it  a  more  concise  form  avoiding  the  explicit  presence  of  6^,  we  express  the 
latter  from  (A. 4).  and  use  the  identity  (28)  and  the  three  lines  following  it: 


^  (P^^’  ^^2-^12,n 


4.n-^12 


Tfius ,  the  desired  term  is 


-2(p^)^2  p'^  4.n-^12 


(A. 7) 


We  are  now  in  a  position  to  write  (A  1.2)  in  a  convenient  form  as  follows: 

<5P  ■  ^,2*-^  ‘  *  '  2  *  *^3  *  ^4  ■  (A- 8) 

where  t  lie  terms  c  ,  c  ,  c  .  and  c  are  given  respectively  by  ( A ,  3 ,  vl ,  6 , 7  )  .  This 

I  fa  O  4 

is  pieeisely  the  final  formulation  (29).  (30a  d).  Although  the  formal  algebraic 
derivation  resulting  in  (A. 1,2)  is  short  and  straightforward  by  comparison  with 
the  geometrir  < onsiderat ions  in  Chapters  2  5.  the  system  (A. 1,2)  is  seen  to  be 

far  from  final.  Its  transformation,  carried  out  with  the  aid  of  relations  from 

Chapters  3  and  5,  is  lengthier  than  the  derivation  of  the  system  (A. 1,2)  itself. 


Tlie  final  formulation  above  has  two  basic  features.  The  first  feature  is 
t)i.it  every  term  on  the  right-hand  side  of  (A. 8)  has  an  easy  geometrical 
interpretation,  and  is  quite  concise  and  tractable.  The  term  c^  depends  on 

,  whereas  .  c^ .  and  c^  all  depend  on  grouped  together  if 

(onvenient  And  the  second  feature  is  that  the  term  c^  has  been  shown  in 

(haptei  C)  to  be  very  significant  and  in  general  mirch  greater  than  the  combined 
term  * '  ,  ‘  '  3  ‘  ^  .j )  •  which,  in  many  applications,  may  be  negligible.  Thus,  a 
sufficiently  accurate  yet  easily  manageable  formulation  appears  to  be 


■^12. n 


(A. 9) 
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